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Abstract 

Experimental evidence has been growing for the existence of both molec- 
ular and exotic dibaryons. The former are dominated by hadron and the 
latter by quark-gluon degrees of freedom. Exotic dihadrons are of partic- 
ular interest because their properties are closely related to the parameters 
and non-perturbative properties of QCD models. Molecular structures pro- 
vide information about confinement. Purely exotic models, such as bag and 
valence quark models, ignore the multi-hadron phase of confinement. Multi- 
configuration quark models and R-matrix boundary condition models better 
represent the medium range repulsive effects leading to confinement, predict- 
ing higher excitation energies of the exotic resonances. The ^5*0 exotic pre- 
dicted at 2.70 GeV mass has been identified in pp spin observables, implying 
that there are several more even and odd parity exotics in the NN system be- 
tween 2.6 and 3.0 GeV. The method also predicts exotic strange and doubly 
strange dibaryons (including a di-lambda unbound by about 0.13 GeV) and 
exotics in baryon and meson channels. The R-matrix formalism also describes 
the properties of molecular structures and enables a classification of these into 
several types. Observed and predicted examples of each type are discussed, 
including the recently observed d' . 
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I. DI-HADRON STRUCTURES 



Just as nuclear bound states and resonances have been important in understanding the 
nature of nuclear forces, the SU(3) structure of the single-hadron {qq and q^) spectrum led 
to the quark model and QCD. It is likely that the multi-hadron spectrum can further reveal 
the properties of QCD as well as of the effective hadron exchange forces. In particular the 
transition from a multi-hadron resonance to single hadrons is also a transition between the 
regions of perturbative asymptotic freedom and non-perturbative confinement in QCD, an 
important process of which we know little. 

But multi-hadron structures have many sources. Some are dominated by long-range 
forces that can be described by hadron exchange potentials (confined quark clusters) and 
others by the short-range quark and gluon degrees of freedom. The former are usually called 
molecular and the latter exotic structures. The long-range dominated structures may arise 
from single- or multi-channel effects. We need to classify the structures by their experimental 
properties in order to deduce the implications for the dynamics. We can learn about these 
relations from a model that adequately includes the physics of both the long- and short- 
range regions. RGM calculations with a sufficiently large number of quark configurations 
approach this behavior and have given useful results for specific cases IQ. The R-matrix 
method [||, which connects the short-range perturbative QCD and the long-range hadronic 
region by a specified boundary condition, has the advantage of a more adequate treatment of 
the non-perturbative region, simpler calculations, and an analytic formulation that enables 
us to classify and describe structures according to the dynamical origin. 

II. QUARK HADRON HYBRID MODEL BY R-MATRIX METHOD 




(a) (b) 
FIG. 1. Two protons at (a) small and (b) large separations. 

Figure |l| represents two components of the pp wave function, (a) being the perturbative 
Dirac quark component with gluon exchange, and component (b) described by a two-proton 
wave function interacting via meson exchange. There are other components like (b) in which 
one or both nucleons are substituted for by the even and odd parity nucleon isobars. Given 
a complete set of internal wave functions ip^ that vanish at Rq, with total energy eigenvalues 
{Wi), R-matrix theory determines the external hadron wave function at energy W, V'vk 



The ^p{r) = rR{r), where R is the radial wave function. 
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by the boundary condition at tq (corresponding to Rq) 



ro^'wiro) = f{W)^y,{ro) with /(W^) = + ^ — ^ (1) 

[the xp^r, f(W), fo, and pi are matrices that couple all isobar components for a given J, /, 
and parity] and 



ft =-ro^U, (2) 



where the are the fractional parentage coefficients of the a channel component of ip-^^ in 
the ipf wave function. Integrating the Schrodinger equation from vq to large r one obtains 
the S-matrix. The positions of the poles and branch points of the S-matrix determine the 
structures (and any bound states) of the two-hadron system. The above model enables us 
to describe the effect of all dynamical components on these singularities, with implications 
for their experimental characteristics. These features are all present in the simple case of 
two coupled channels with different thresholds (such as NN ^ AA), with the neglect of 
the off-diagonal potential. Then the scattering of an incoming channel of orbital angular 
momentum L is described by H 



where 



Where k and k' are the relativistic relative momenta in the L and L' channels, 

and the J^^'^ are the outgoing and incoming Jost functions, approaching the spherical 
Hankel or the Coulomb functions at large range. Their analytic properties determine that 
the 6l are real for imaginary k (below threshold) and are complex for real k. For a negligible 
strength potential 



6'L(K,ro) = — TjT- ^ >L + tkro. (6) 

hl'{kro) 



III. SOURCES OF SCATTERING AMPLITUDE STRUCTURES 

Poles of the S-matrix in the complex energy {W) plane will produce significant structure 
in the energy dependence of scattering amplitudes if they are close enough to the physical 
cut. Or, if they are on the real W axis, then they must be below the minimum strongly 
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(e) 



FIG. 2. Lowest order diagrams representing the various mechanisms that can cause structure 
(resonances or bound states) in two-hadron interactions. For definiteness, baryons ( — ) are shown 
interacting via mesons ( — ) or a core (o) representing the average effect of short-range meson and 
gluon exchange, and intermediate states with baryon isobars (=) or quarks ( — ) and gluons('wv). 

coupled two-particle mass [k imaginary for all channels) and represent a bound state. Also 
significant structure may be caused by inelastic branch points [k' = for a coupled channel) 
with a sufficiently strong discontinuity, representing strong coupling between channels. 

The sources and character of each type of structure can be inferred from Eqs. A 
pole of the S-matrix requires a zero of DEN(Vl^) = feff(W) + 6'l(A;, ro), and an observable 
structure in the continuum further requires that the W dependence of DEN be sufficiently 
rapid in the vicinity of the zero. The diagonal potential in channel 1 is represented by 
ro). Eq. I shows that if fesiyV) < —L (see Fig.^(a)) there will be a bound state in 
the absence of a potential. The attraction is then due either to the core fii or to the 
coupled channels in /eg. If the potential is attractive (see Fig.||(b)), Ol~ L will be negative 
below threshold and there may be a bound state for /eg > —L or even for a repulsive core 
/cff > L + 1. From Eq. ^ one can infer that in the continuum (real fc), the pole will not be near 
the physical cut (producing a narrow resonance) unless k is small (an "anti-resonance" ) or 
the potential is nonmonotonic producing a potential barrier and rapid variation of ^^(fc, tq). 
Alternatively the rapid variation may be in the /eg term, which by Eq.|^ requires either 
that the coupling be large (a large potential coupling will have the same effect - see 

Fig.H(c)), or that the denominator DEN' = + ^L'{k',rQ) vanish (see Fig.||(d)) for 

k'^ < 0. The former causes a resonance below the second channel threshold or even a bound 
state due entirely to interchannel coupling §]. For small the associated pole is 

far from the physical cut and the structure is an L-th order cusp produced by the inelastic 
branch point. As increases the pole comes nearer, and the structure forms into a 

broad resonance below the threshold going further down in energy and becoming narrower 
with further increase of The latter case (vanishing DEN') is a state that would be 

bound if there were no coupling to the first channel, but leaks into that channel through the 



coupling. The smaller the coupling the narrower the resultant "Dalitz-Tuan" resonance . 
The position of the resonance stays near the zero of DEN'. If the channels are only coupled 
by weak or electromagnetic forces, the result will be a very narrow resonance. None of the 
above structures are importantly affected by the positions or residues of the poles due to 
the interior quark states. 

However, in the vicinity of quark configuration energies (see Figs.[l| and 0(e)) the /- 
matrix poles cause a rapid variation of fcs- The positions of these resonances are near the 
quark state energies, and the splittings for a given quark configuration are determined by 
the one-gluon exchange matrix elements. The width of the resultant structure depends on 
(Eq.[^), which is decreased from the scale of r^^ by the fractional parentage coefficients. 
This results in widths of ~ 50 MeV to ~ 100 MeV, easily detected with the resolution 
of multi-GeV beams. It is important to note that any very narrow resonances observed 
experimentally are more readily described by molecular "Dalitz-Tuan" resonances than by 
exotic resonances induced by quark degrees of freedom. 

IV. HADRONIC EXAMPLES OF VARIOUS STRUCTURES 

There are many examples of potential bound molecular states, the deuteron and other 
nuclei among them. Barrier potential resonance examples are represented by a unstable and 
fissioning nuclei. The and ^F^ NN states near the A^A threshold are well understood 
as driven by the A^A threshold branch point. The yg*(1405) and the S* are caused 
by the coupling to K~p and KK systems, respectively; either through the strong coupling 
or by the "Dalitz-Tuan" formation of a quasi-bound state |Q. The recently observed |^ 
(i'(2060) has a width < 1 MeV, and could be described as a bound state of a nucleon 
and an odd-parity nuclear isobar I{J^) = 0{0~) coupled electromagnetically to an NNtt 
state. 

The spectra of the quark exotic [g(lsi)]^, [g(lsi )]^g(lpi ), [g(lsi )]^s(lsi ), 
[g(lsi )]''[s(lsi )]^, and [g(lsi )]'^s(lsi ) quark configurations have been predicted by the R- 
matrix method @J^,|] {q represents u or d quarks). Constituent quark calculations predict 
similar spectra, at about 100 MeV less mass. The lowest mass 1 = 1 NN state is predicted 
at 2.7 GeV/c^ mass. Several experiments have observed structure at this energy, the most 
significant of which (see Fig. ^ is the energy dependence of the spin correlation Ann{pp) [Q. 
The prediction made in 1987 is in good agreement with the data. Some recent np data 
is consistent with the lowest predicted / = state at 2.63 MeV 0], but data is required in 
finer energy steps. We await further experiments in the NN, AN — SA^, AA — SN, and KN 
systems. 
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Aoonn(90) vs. Energy (May/June 1994 running period only) 
(Three 2-degree angle bins around 90 degrees were averaged at each energy) 
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Theory of Ref. [2] 



No background subtraction was done 
before producing these points. 
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